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ABSTRACT 

Evaluation  of  a  real  function  f  on  an  interval  X  using  interval  arithmetic 
yields  an  interval  extension  F(X)  containing  the  range  R(f;X)  of  f  on  X. 
Unfortunately,  F(X)  is  sometimes  excessively  wider  than  R(f;X).  Evaluation  of 
f  t  C1  by  interval  differentiation  arithmetic  gives  F(X)  and  the  extension  F'(X) 
of  f'  on  X.  If  F'(X)  >0  (or  F*(X)  <  0),  then  f  is  monotone  on  X  *  [a,b] ,  and 
R(f;X)  =  [f(a),f(b)]  (or  R(f;X)  =  [f (b) , f (a ) ] ) ,  giving  improved  bounds  for 
R(f;X).  If  0  e  int(F'(X)),  then  X  is  divided  into  subintervals  on  which  f  is 
either  guaranteed  to  be  monotone  or  has  possible  extremal  points. 

A  Pascal-SC  program  for  this  single  algorithm  is  given,  and  numerical 
results  are  presented.  As  a  byproduct  of  the  computation,  possible  extremal 
points  of  f  are  isolated.  \ 


AMS  (MOS)  Subject  Classifications:  65G10,  65K10 

Key  words:  Range  of  Functions/  Extremal  Points,'  Automatic  Differentiation," 
Interval  Computation/  Optimization,  Pascal-SC 
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SIGNIFICANCE  AND  EXPLANATION 


Interval  computation  provides  ways  to  obtain  inclusions  for  the  range  of  a 
function  on  an  interval  automatically,  taking  into  account  the  effects  of 
roundoff  error  in  actual  computer  evaluation  of  the  function*  However,  simply 
"plugging  in  an  interval"  for  the  independent  variable  and  using  interval 
arithmetic  to  evaluate  the  function  often  leads  to  inclusions  of  the  range  of  a 
function  which  are  too  large  by  orders  of  magnitude,  and  useless  for  practical 
purposes.  Because  of  this  failure  of  naive  application  of  interval  methods,  a 
number  of  efficient  methods  for  calculating  an  inclusion  of  the  range  of  a 
function  have  been  developed,  most  of  which  are  based  on  the  use  of  a  mean-value 
or  centered  form  for  the  interval  extension  of  a  function.  The  algorithm 
presented  in  this  paper  takes  a  different  approach,  based  on  the  fact  that  the 
range  of  a  monotone  function  on  an  interval  can  be  computed  by  evaluating  the 
function  at  the  endpoints  of  the  interval  (with  directed  rounding  to  insure 
inclusion).  Interval  differentiation  arithmetic  is  used  to  evaluate  an  interval 
inclusion  F'(X)  of  the  range  of  the  derivative  of  the  function  on  the  interval 
X.  Even  though  this  inclusion  may  be  crude,  the  function  is  guaranteed  to  be 
monotone  on  X  if  0  is  not  contained  in  the  interior  of  F ' (X) ,  and  hence  its 
range  can  be  calculated  accurately.  If  0  is  in  the  interior  of  F'(X),  then  X  is 
subdivided  into  subintervals,  on  which  the  function  is  either  monotone  or 
possibly  has  an  extremal  point.  By  taking  the  interval  hull  of  the  ranges  of 
the  function  on  all  subintervals,  an  improved  inclusion  of  the  range  of  the 
function  is  obtained.  As  a  byproduct  of  the  computation,  possible  extremal 
points  of  the  function  are  isolated  in  the  subintervals  on  which  0  is  in  the 
interior  of  the  interval  extension  of  the  derivative  of  the  function.  If  there 
are  no  such  intervals,  then  the  function  is  piecewise  monotone,  and  an  accurate 
inclusion  of  its  range  is  obtained  by  the  algorithm. 


A  Pascal-SC  program  is  given  which  implements  the  algorithm,  and  numerical 
results  are  presented.  The  text  of  this  report  was  presented  as  an  invited 
address  at  the  International  Interval  Symposium  1985,  held  at  the  University  of 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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1.  Ranges  of  functions.  The  range  of  a  real  function  fsD  C*  -*•  R  on  a  set  X  C  D  Is 


R(fiX)  -  {f(x)  |  x  e  X}. 


In  case  X  «  [a,b]  is  a  closed,  bounded  interval  and  f  is  continuous,  then  R(f;X)  will  also 
be  an  interval  of  the  same  kind.  Closed,  bounded  intervals  will  be  referred  to  simply  as 
intervals,  and  the  set  of  such  intervals  will  be  denoted  by  IR. 


ly.'.vV", 


A  fundamental  problem  of  interval  analysis  is  the  calculation  of  R(f;X)  or  at  least  a 
good  approximation  to  it.  If  f  is  defined  in  terms  of  arithmetic  operations  and  functions 
with  known  interval  extensions,  then  straightforward  use  of  interval  computation  gives  an 
interval  extension  F  of  f  such  that 


R ( f ; X )  C  F(X) 


■•-Wa'J 


for  X  C  D.  This  calculation  has  the  advantage  of  being  completely  automatic,  and  does  not 
require  knowledge  of  special  properties  of  f.  Unfortunately,  F(X)  can  be  such  a  gross 
overestimation  of  R(f|X)  in  certain  cases  that  it  is  useless  for  practical  purposes. 
Furthermore,  the  quality  of  F(X)  as  an  approximation  to  R(f;X)  is  generally  unknown. 
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A  number  of  methods  have  been  developed  for  obtaining  better  approximations  to  R(fiX), 
starting  with  the  work  of  Moore  [1].  The  recent  book  by  Ratschek  and  Rokne  [51  describes  a 
number  of  these  techniques,  and  gives  a  substantial  bibliography.  Most  of  the  apprc  iches 
to  this  problem  are  based  on  transformation  of  F,  usually  into  centered  or  mean-valu  orms 
[11,  [51.  The  method  given  in  this  paper  applied  to  continuously  differentiable  tions 
f,  and  makes  use  of  information  about  the  monotonicity  of  f  obtained  by  the  p.  i  of 
automatic  differentiation  [2] . 

2.  Monotone  functions.  If  the  function  f  is  nondecreasing  on  X,  then  R(f  )  is 
simply 

R(fiX)  -  [f (a) ,  f(b)].  (2.1) 

Similarly,  if  f  is  nonincreasing  on  X,  then 

R(f;X)  -  (f(b),  f(a) ) .  (2.2) 

Thus,  the  range  of  monotone  functions  can  be  determined  by  calculating  only  two  function 
values.  In  actual  practice,  of  course,  downward  rounding  of  the  lower  endpoint  and  upward 
rounding  of  the  upper  endpoint  gives  an  interval  inclusion  of  R(f;X)  which  is  slightly 
wider  than  the  exact  range.  For  the  time  being,  it  will  be  assumed  that  function  values 
are  computed  exactly. 

A  sufficient  condition  for  (2.1)  to  hold  for  differentiable  f  is  that 

f'(x)  >0,  a  <  x  <  b,  (2.3) 

and  similarly  (2.2)  holds  if  f'(x)  <  0  on  X.  Furthermore,  suppose  that  f  is  continuously 
differentiable,  and  F'  denotes  an  interval  extension  of  f'  obtained  by  interval 
computation.  If  F'(X)  >  0  (F'(X)  <  0),  it  follows  that  f  is  nondecreasing  (nonincreasing) 
on  X,  and  R(f;X)  can  be  calculated  directly  by  (2.1)  or  (2.2),  respectively. 


The  additional  information  about  the  derivative  of  f  needed  above  can  also  be  obtained 
automatically.  The  values  of  F(X)  and  F'{X)  can  be  computed  by  using  interval 
differentiation  arithmetic,  as  described  below.  All  that  is  required  is  a  formula  or 
subroutine  for  f;  no  symbolic  differentiation  is  necessary.  If  necessary,  a  bisection 
procedure  can  be  applied  to  the  interval  X  to  find  subintervals  on  which  f  can  be 
guaranteed  to  be  monotone.  The  resulting  algorithm  provides  either  the  exact  value  of 
R(f;X),  or  else  an  inclusion  of  R(fiX)  which  is  better  in  general  than  F(X). 

3. _ Real  differentiation  arithmetic.  It  is  convenient  to  define  interval 

differentiation  arithmetic  as  an  extension  of  real  differentiation  arithmetic.  This 
arithmetic  can  be  used  to  calculate  the  values  of  functions  and  their  derivatives 
automatically,  without  symbolics  or  numerical  approximations  [4).  Like  interval 
arithmetic,  real  differentiation  arithmetic  is  an  ordered-pair  arithmetic,  with  elements  U 
=  (u,u'),  V  =  (v,v‘),  ...  e  R .  The  rules  for  this  arithmetic  are: 


U  +  V  -  (u,u')  +  (v,v')  »  (u  +  v,  u'+  v'),  (3.1) 

0  -  V  »  ( u , u 1 )  -  (v,v')  «  (u  -  v,  u*  -  v'),  (3.2) 

U«v  »  (u,u' )  *(v,v' )  «  (u«v,  u»v'  +  v*u'),  (3.3) 

U/V  »  ( u , u ' )/(v,v' )  =  (u/v,  (u'  -  (u/v) • v 1 )/v),  v  +  0.  (3.4) 


The  arithmetic  defined  in  this  way  forms  a  division  ring  with  identity,  and  will  be  denoted 
by  D.  If  the  first  element  of  each  operand  pair  is  interpreted  as  a  function  value,  and 
the  second  as  a  derivative  value,  then  the  first  element  of  the  result  corresponds  to  the 
evaluation  of  the  operation,  and  the  second  to  the  evaluation  of  its  derivative,  according 
to  the  well-known  rules  of  calculus.  If  real  numbers  c  are  identified  with  the  pair3 
(c,0),  then  it  follows  from  the  chain  rule  of  calculus  that 

f ( ( x, 1 ) )  =  (f(x),  f'(x)),  (3.5) 


that  is,  the  rules  of  differentiation  arithmetic  will  automatically  give  both  the  value  and 
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the  value  of  the  derivative  of  a  rational  function  f.  More  generally/  the  chain  rule  gives 


fUu/U1))  =  (f(u),  u'.f{u)). 


(3.6) 


which  allows  the  definition  of  standard  functions  in  D,  for  example. 


> 

■  „U  _  ( u ,  u  * ) 

e  =  e 

In  U  =  ln(u,u' ) 


,  U  ,  U\ 

(e  ,  u'«e  ), 
(In  u,  u'/u). 


(3.7) 

(3.8) 


I 

and  so  on.  The  combination  of  arithmetic  operations  and  standard  functions  will  be  called 
a  computational  system  for  differentiation  arithmetic.  It  is  simple  to  program  such  a 
computational  system,  particulary  in  a  language  such  as  Pascal-SC,  which  permits  definition 
of  operators  and  functions  for  various  data  types  [3] . 

4.  Interval  differentiation  arithmetic.  Interval  differentiation  arithmetic  is 
defined  by  the  same  rules  as  real  differentiation  arithmetic,  starting  with  pairs  of 
intervals  instead  of  real  numbers,  and  using  interval  arithmetic  instead  of  real  arithmetic 
inside  the  parentheses  on  the  right  sides  of  (3.1 )-(3.3) .  With  interval  extensions  of 
standard  functions,  the  definitions  (3.7),  (3.8)  and  so  on  are  used  to  construct  a 
computational  system  for  interval  differentiation  arithmetic.  Once  again,  such  a  system  is 
easy  to  program  in  Pascal-SC,  which  supports  interval  arithmetic  as  well  as  operator  and 
function  definitions  for  various  data  types  [3]. 

The  analog  to  (3.5)  in  interval  differentiation  arithmetic  is 


F( (X,  (1, 1) ))  =  { F (X ) ,  F' (X)  ).  (4.1) 

Thus,  by  a  direct  evaluation  process  in  this  arithmetic,  interval  inclusions  F(X)  of  R(f;X) 
and  F '  ( X )  of  R(f';X)  can  both  be  obtained  automatically.  Here,  even  if  F'(X)  is  a  crude 
approximation  to  the  range  of  f*  on  X,  the  conditions  F'(X)  >  0  or  F’(X)  <  0  are  sufficient 
to  guarantee  the  monotonicity  of  f,  and  if  f  is  monotone,  then  its  range  can  be  calculated 


exactly  by  (2.1)  or  (2.2).  This  observation  is  the  basis  of  the  algorithm  described  in  the 
next  section. 

5.  to  algorithm  for  range  calculation.  Of  course,  if  the  calculation  of  F'(X)  shows 
that  f  is  monotone  on  the  entire  interval  X,  then  R(f ;X)  can  be  calculated  at  once. 
Otherwise,  X  will  be  partitioned  into  subinterval,  and  either  R(f)X)  or  an  approximation  to 
it  will  be  constructed.  Let  a  given  list  of  n  subintervals  of  X  be  denoted  by  Lr  «  {X^, 
X_,...,  XR),  and  suppose  that  R  C  R(ftX)  is  known.  On  each  subinterval  X^,  either  F(X^)  C 
R,  in  which  case  R(f;X^)  makes  no  additional  contribution  to  R(f|X),  or  f  is  monotone,  in 
which  case  its  range  can  be  computed  directly  and  R  updated,  or  else  0  is  an  interior  point 
of  F’(X^),  in  which  case  may  contain  a  critical  point  of  f.  In  the  latter  case,  X^  can 
be  bisected  and  the  resulting  subintervals  put  on  a  now  list  for  further  examination.  In 
order  for  the  algorithm  to  terminate  in  a  finite  number  of  steps,  a  lower  bound  $  is  put  on 
the  widths  of  the  subintervals  to  be  considered,  and  an  upper  bound  N  is  placed  on  the 
number  of  subintervals  to  be  saved  for  further  examination.  For  convenience,  if  Y,Z  are 
intervals,  then  Y  ++  Z  will  denote  the  interval  hull  of  Y  and  Z,  that  is,  the  smallest 
interval  which  contains  both  Y  and  Z. 

The  algorithm  consists  of  the  following  stepsi 

1*.  (Initialisation)  Taka  X1  :«  X,  L,  i«  {X,},  R  (f(x),f(x)],  where  x  is  some 
point  in  X. 

2* •  (Iteration)  For  1  «  1,...,n,  compute  (F(X^),  F'lX^)) 

(a)  If  FtX^  C  R,  then  discard  X^ 

(b)  If  F 1 ( X^ )  >  0  or  F'tXj)  <  0,  then  compute  R  :»  R  ++  Rtf.X^)  and  discard  X^ 

(c)  Otherwise,  retain  . 

3°.  (Termination  or  continuation)  Denote  the  list  of  retained  intervals  by  Lf. 

(a)  If  Lf  is  empty,  then  the  algorithm  terminates  with  the  exact  value 

R  -  R(f IX)  (5.  1 ) 


of  the  range  of  f  on  X 


(b)  If  r  >  N  or  w(x, )  <  z,  then  the  algorithm  terminates  with  the  overestimate 

R  R  ++  F(X j )  ++...++  F(Xr)  D  R(f)X)  (5.2) 

of  the  range  of  f  on  X. 

(c)  Otherwise,  each  subinterval  in  Lr  is  bisected  to  form  a  new  list  with  n  » 
2r,  and  the  algorithm  returns  to  step  2*. 

6.  Remarks.  The  algorithm  given  in  the  previous  section  will  terminate  in  a  finite 
number  of  steps  with  either  the  exact  value  of  R(f;X)  or  an  overestimate  which  is  never 
worse  than 

R  -  FfX^  ++  ...  ++  F(Xn)  3  R(f)X).  (6.1) 

In  general,  (6.1)  is  a  better  approximation  to  R(f,-X)  than  F(X)  because  of  the  convergence 
of  united  extensions  to  the  range  of  a  continuous  function  [1]. 

As  a  byproduct  of  the  calculation  when  an  overestimate  is  produced,  the  intervals 
X1,...,Xr  which  are  retained  at  the  final  step  may  contain  critical  points  of  f,  that  is, 
points  at  which  f'(x)  -  0.  This  information  may  be  useful  in  optimization  problems. 
Furthermore,  if  the  list  of  retained  intervals  is  nonempty,  then  the  value  R  3  R(f;X) 
returned  by  the  algorithm  is  definitely  known  to  be  an  overestimate,  while  if  the  list  of 
retained  algorithms  is  empty,  then  this  value  is  exact  (modulo  outward  rounding).  Thus, 
the  algorithm  itself  indicates  the  type  of  result  (exact  or  an  overestimate)  it  obtains. 
The  knowledge  that  R  is  an  overestimate  and  the  list  of  retained  intervals  can  be  used  to 
refine  the  calculation  of  R(f;X)  further,  if  desired.  Some  idea  of  the  quality  of  the 
overestimate  can  be  obtained  by  comparing  the  value  of  R  before  calculating  (5.2)  with  the 
final  result. 

7.  Numerical  results.  Numerical  results  were  computed  for  the  following  functions, 
using  the  Pascal-SC  program  given  in  the  following  section. 


f^X)  »  X  -  X 


(7.1) 


f2(x)  »  x»x. 


f  j(x) 


(x  -  1)-(x  »  3) 
(x  +  2) 


f4(x)  »  x/x. 


(7.2) 

(7.3) 

(7.4) 


a.  For  X  «  [a,b],  the  naive  interval  extension  F^tX)  -  X  -  X  of  £1  gives  F^ta.b])  - 
(a  -  b,  b  -  a] ,  while  the  algorithm  gives  R  «  (0,0)  -  R(f.,iX)  for  arbitrary  X. 

b.  For  symmetric  intervals  X  «  [-s,s),  the  algorithm  gives  the  exact  value  R  «  (0,s2) 

“  R(f2;  [-8,8]),  while  F2([-s,s))  »  (-s,sl • (-s,s)  *  t-s2,s2].  In  case  X  »  (-r,s)  is 

nonsymmetric  interval  containing  0,  the  result  of  the  algorithm  can  be  of  the  form  R  « 

(-E,  max{r  ,s  },  where  e  >  0  is  small,  with  a  message  that  a  small  interval  containing  0 
can  contain  a  critical  point  of  fj.  For  example,  for  X  -  (-7,8),  one  has 

F2(X)  -  X-X  -  (-56,  64)  (7.5) 


while  the  algorithm  gives 


R  -  (-3.1x10_,e,  641  (7.6) 

with  a  notation  that  there  may  be  a  critical  point  of  f2  in  the  retained  interval 

_>g 

[-1.63x10  ,  1.87x10  '] .  In  all  other  cases,  the  algorithm  gives  the  exact  result.  Even 

if  X  is  nonsyirmetric  about  0,  the  algorithm  will  give  the  correct  result  if  0  is  a 
bisection  point. 

c.  The  function  f^  is  actually  monotone  increasing,  but  has  a  pole  at 
x  =  -2.  The  algorithm  will  sense  the  monotonicity  of  f^  and  give  correct  results  if  X  is 
subdivided  a  sufficient  number  of  times.  The  results  are  such  better  than  the  naive 
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interval  extension  F 3<X)  •  (X  -  1 )  * (X  +  3)/(X  +  2)  when  one  of  the  endpoints  of  X  is  close 
to  -2.  For  example,  for  X  «  [-1.9,  98) , 


F3(X)  -  [-2929,  97970], 


while  the  algorithm  gives 


R  -  [-31.9,  97.97], 


For  X  -  [-1.999999,  98]  , 


while  the  algorithm  gives 


F,(X)  -  [-3.03x10°,  9.797x10  ], 


R  -  [-3000002,  97.97). 


Finally,  for  X  »  [-1.99999999999,  98],  which  has  a  lower  endpoint  as  close  to  -2  as 

possible  in  12-digit  decimal  arithmetic,  one  gets 


Fj(X)  =  [-3.03x1013,  9.797x1014], 


•  •] 


while  the  algorithm  gives 


R  »  [-300000000002,  97.97) 


d.  The  algorithm  does  not  give  good  results  for  f4(x)  “  x/x,  because  it  determines 


that  every  subinterval  of  X  possibly  contains  a  critical  point  of  f^  (which  in  fact  is 
true,  since  f4(x)  =  1  is  constant,  and  f4‘(x)  =  0).  Thus,  the  algorithm  computes  R  *  [1,1] 
initially,  and  the  final  value  is  determined  only  by  the  united  extension  (6.1).  of 
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course,  the  result  is  generally  better  than  the  naive  interval  extension  F4(X)  «  X/X 

evaluated  on  the  entire  interval  X,  but  is  still  usually  a  gross  overestimate.  For 
example,  for  X  «  [0.002,  2]. 

F4 (X )  -  [0.001,  1000]  (7.13) 


while  the  algorithm  gives 


R  -  [0.203,  4.903],  (7.14) 

which  is  still  not  a  very  good  approximation  to  [1,1],  even  though  it  is  much  better  than 
(7.13).  Of  course,  the  user  is  warned  that  the  result  may  not  be  good  by  the  fact  that  all 
subintervals  are  retained.  Other  methods  usually  give  no  warning  when  gross  overestimates 
are  produced.  One  way  to  improve  the  algorithm  in  this  case,  since  interval  extensions  of 
derivatives  are  available,  would  be  to  use  mean-value  forms 

F(Xt)  -  m(XJL)  +  F'(XjL).(X  -  m(Xi))  (7.15) 

to  obtain  interval  extensions  F  of  f  on  subintervals  X^,  instead  of  obtaining  them  by 
straightforward  evaluation. 

8.  A  Pascal-SC  program.  The  program  written  below  was  designed  to  be  general,  so 
that  the  user  needs  to  supply  only  subroutines  for  evaluation  of  the  function  f  in  ordinary 
interval  arithmetic  (IFEVAL)  and  in  interval  differentiation  arithmetic  (IDFEVAL) .  The 
source  code  for  these  subroutines  should  be  located  in  the  files  FEVAL.FUN.  Examples  of 
these  subroutines  for  the  functions  discussed  in  §7  are  given  in  §10. 

The  operators  for  interval  differentiation  arithmetic  given  in  §9  include  only  the 
basic  arithmetic  operators  for  type  IDERIV.  For  a  complete  computational  system,  operators 
for  mixed  arithmetic  between  types  INTEGER,  REAL,  and  IDERIV  should  be  included,  as  well  as 


standard  functions  [3]. 


V.TP.TJf.  r.  *  ■  "  I  ^  •  '  - 


The  number  of  subintervals  allowed  in  a  list  is  set  by  the  constant  DIM  in  the 
program,  which  can  be  changed  by  the  user.  The  size  of  the  smallest  subintervals  is 
similarly  controlled  by  the  number  LIMIT  of  bisections  allowed.  Thus,  if  LIMIT  “  L,  then 

&  ■=  2"L.w(X),  (8.1) 

where  w(X)  =  b  -  a  is  the  width  of  the  original  interval  X  *  [a,b]  ,  The  source  code  for 
the  Pascal-SC  program  follows: 


PROGRAM  IRANGE ( INPUT, OUTPUT ) ; 

CONST  DIM  =  256j  (*  Maximum  number  of  subintervals  *) 

LIMIT  =  32;  (*  Maximum  number  of  bisections  *) 

TYPE  INTERVAL  =  RECORD  INF, SUP  :  REAL  END  ; 

IDERIV  -  RECORD  X, PRIME:  INTERVAL  END: 

DIMTYPE  =  1..DIM; 

STACKTYPE  «  RECORD  INT : INTERVAL;FUN : IDERIV  END: 

VAR  X,RF, BEST, WORST:  INTERVAL; 

F:  IDERIV; 

I,NA.NB,LIM:  INTEGER; 

A,B:  ARRAY [DIMTYPE]OF  STACKTYPE;  (*  A  is  the  list  of  intervals  to 

be  examined,  B  is  the  list  of 
retained  intervals  *) 

MX:  REAL; 

^INCLUDE  INTERVAL. PAK;  (*  Makes  interval  arithmetic  available  ») 

$INCLUDE  IDERV.PAK;  (*  Interval  differentiation  arithmetic  *) 

PROCEDURE  IOUT(X:  INTERVAL);  (*  Prints  endpoints  in  standard  format  *) 
BEGIN 

WRITE ( ’ [ 1 ,X. INF, • , • ,X.SUP, • ] ' ); 

END; 
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SXNCLUDE  FBVAL.FUN; 


(*  Evaluation  of  the  function  in  interval  and 


and  interval  differentiation  arithmetic  *) 
FUNCTION  RMF(L,G:  REAL):  INTERVAL: 

(*  Bounds  the  range  of  a  monotone  function  which  assumes  its  least 
value  at  L  and  its  greatest  value  at  G.  *) 

VAR  D,U:  INTERVAL; 

BEGIN 

D:-INTPT(L);U:-INTPT(G)j 
D:«IFEVAL(D) ;U:-IFEVAL(U)| 

D.SUP:«U.SUP; 

RMF:  *D 
END; 

function  MID(X:  INTERVAL):  REAL;  <*  Calculates  midpoint  of  an  interval  *) 
VAR  A,B:  ARRA5f(T.  .2JOF  REAL; 

BEGIN 

A[1] :»X.INF;B[1] :-0.5; 

A [2] :«X.SUP;B[2J :-0.5; 

MID:«SCALP(A,B,0) 

END; 

BEGIN  (*  Program  IRANGE  *) 

WRITELN( 'Enter  initial  interval  X:'); 

IR£AD( INPUT* X) ; 

WRITE ('  X-  '  );IOUT(X);WRITELN; 

F:«IDFEVAL(X); 

WORST :-F.X; 

IF  (F. PRIME. INF  >«  0)  THEN  RF:»RMF(X. INF, X. SUP) 

ELSE  IF  (F. PRIME. SUP  <«  0)  THEN  RF: “RMF (X. SUP,X. INF ) 


BEGIN  I*  Mi  not  fflonotona  * ) 

NA:-1|LIM:-0) 

AC11.INTi-X|A[1] .FUNi-F; 

MX:-MID(X); 

X: “INTPT(MX ) I 
BEST : "IFEVAL (X ) ) 

WHILE  ((NA  >  0)  AND  (NA  <-  DIM  DIV  2)  AND  (LIM  <  LIMIT))  DO 
BEGIN  (*  WHILE  *) 

LIM : "LIM+ 1 ;  NB : «0 1 
FOR  I :  — 1  TO  NA  DO 
BEGIN  (*  STACK  B  *) 

MX:-MID(A[I]  .INT)| 

NB:«NB+1; 

B[NB]  .INT.INF:-A[I] .INT.INF) 

B[NB] .INT.SUP:-MX) 

B[NB]  .FUN:«IDFEVAL<B[NB]  .INT); 

NBi-NB+1; 

B [NB] . INT* INF: “MX | 

B  [NBj . INT. SUP: “A [I] .INT. SUP) 

B  {NB]  .  FUN:**IDFEVAL(B  [NB]  •  INT )  | 

END;  (*  STACK  B  *) 

NA:«0; 

FOR  I : “1  TO  NB  DO 
BEGIN  (*  UNSTACK  B  *) 

IF  NOT  ( B [ I ] .FUN.X  <-  BEST) 

THEN  IF  (B [I] .FUN. PRIME. INF  >-  0) 

THEN  BEST : -BEST+*RMF ( B (I ) . INT . INF , B [I ] . INT . SUP ) 

ELSE  IF  (B [I] .FUN. PRIME. SUP  <-  0) 

THEN  BEST:«BEST+*RMF(B(I] . INT. SUP, B [I] .INT. INF) 


-12- 


BEGIN  (*  RESTACK  A  •) 


NA:-NA+1)A[NA]  :-B{I] 

END)  (*  RESTACX  A  *> 

END)  (*  UNSTACK  B  *) 

RF :  “BEST ) 

FOR  1 1  —  I  TO  NA  DO  RF:-RF+*A[I)  . FUN.Xl 
END)  (*  WHILE  *) 

IF  NA  >  0  THEN 
BEGIN  (*  NA  >  0  •) 

RF : -BEST) 

WRITELN( 'Function  may  have  critical  points  ini')) 
FOR  I :  "1  TO  NA  DO 


BEGIN 

WRITE { 'A -  • ))IOUT(A[I] . INT)|WRITELN) 
RF«-RF+«A[I] -FUN.X 


END) 

END)  <•  NA  >  0  •) 

END)  (*  F  ia  not  monotone  *) 

WRITELNl 'Naive  interval  arithmetic  gives:')) 
WRITE ( '  F(X)  -  '  )  )  IOUT ( WORST )  )WRITEUI ) 
WRITELN ( ’The  algorithm  gives:')) 

WRITE (  '  F(X)  -  '  )  )  IOUT ( RF ) ) WRITELN 
END.  (*  Program  I RANGE  *) 


eratora  for  interval  differentiation  arithmetic.  The  six  basic  unary  and 


binary  arithmetic  operators  for  type  IDERIV  are  located  in  the  file  IDERIV.PAK,  which  also 


includes  the  call  to  the  interval  library  for  the  function  ISCALP  to  confute  the  interval 


scalar  product 


TYPE  I  VECTOR  -  ARRAY!  1 ..  2JOP  INTERVAL) 

FUNCTION  I  SC  ALP  (  VAR  A,B:  IVECTORl  DIM)  INTEGER):  INTERVAL) 
EXTERNAL  88;  (*  Interval  scalar  product  * ) 

OPERATOR  +  (U:  IDERIV)  RES:  IDERIV; 

BEGIN 
RES :  *U 
END; 

OPERATOR  -  (U:  IDERIV)  RES:  IDERIV) 

BEGIN 
U.X:— U.X) 

U .  PRIME :  —  U .  PRIME ; 

RES :  “U 
END; 

OPERATOR  +  (U.V:  IDERIV)  RES:  IDERIV) 

BEGIN 

U.X:«U.X+V.X) 

U .  PRIME :  «U.  PRIME +V.  PRIME; 

RES:-U 

END; 


OPERATOR  -  (U,V:  IDERIV)  RES:  IDERIV; 
BEGIN 

U.X:”U.X-V.X; 

U . PRIME : ”U . PRIME -V . PR IME ) 

RES:  -U 


aS&jgtti 


OPERATOR  •  (U,V:  IDERIV)  RES:  IDERIV; 


VAR  A,Bs  I  VECTOR) 


A[1]  :-U.X;B[1]  :-V. PRIME: 
A  [2]  :«V.X;B[2J  :-V. PRIME: 


U .  PRIME :  -I  SCALP  (  A ,  B ,  2  )  : 


U.X:“U.X*V.X; 


OPERATOR  /  (O, V:  IDERIV)  RES:  IDERIV; 


ti.tifk 


•‘-VW 


•  i 


VAR  A,B:  I  VECTOR: 


C:  IDERIV; 


.  tJ 


C.Xs-U.X/V.X; 


A [ 1 ] :-INTPT(1):B[1] :-U. PRIME; 
A[2) C.X;B[2] :-V. PRIME: 


C.PRIME:-ISCALP(A,B,2)/V.X; 


•\  •  .  - 
V.  t 


IFEVAL  :»  X  -  X 


A- *  * 


to.  Example  function  subroutines.  (Content*  of  the  file  FEVAL.FUN.) 


(a)  f^x)  -  x  -  x. 


FUNCTION  IFEVAL (X:  INTERVAL):  INTERVAL; 


’.VV.N 


-i.  . ■  -k.  i 


FUNCTION  IDFEVAL  (X:  IDERIV):  IDERIV: 

BEGIN 

IDFBVAL  X  -  X 
END; 

(b)  f 2 ( x )  -  x*x. 

FUNCTION  IFEVAL (X:  INTERVAL ) :  INTERVAL; 

BEGIN 

IFEVAL  :»  X*X 
END; 

FUNCTION  IDFBVAL (Xs  IDERIV):  IDERIV; 

BEGIN 

IDFEVAL  X*X 
END 

(c)  f3(x)  -  (x  -  1 )  - (x  +  3)/(x  +  2). 

FUNCTION  IFEVAL (X!  INTERVAL);  INTERVAL; 

BEGIN 

IFEVAL  (X  -  1 ) * (X  +  3)/(X  +  2) 

END; 

FUNCTION  IDFEVAL (X:  IDERIV):  IDERIV; 

VAR  ONE, TWO. THREE:  IDERIV; 

BEGIN 

ONE.X  :»  INTPT ( 1 ) ; ONE .  PRIME  :-  INTPT(O); 
TWO.X  :-  INTPT(2) ;TWO. PRIME  :«  INTPT(O); 
THREE. X  :-  INTPT  ( 3  );  THREE .  PRIM  :-  INTPT(O); 
IDFEVAL  :-  (X  -  ONE)»(X  +  THREE )/(X  +  TWO) 
END; 


<d)  f  (x)  *  x/x 


FUNCTION  inVU(Xl  INTERVAL ) !  INTERVAL) 


I 


t 


BEGIN 

IFEVAL  :»  X/X 
END) 

FUNCTION  IDFEVALtX:  IDERIV) !  IDERIV) 

BEGIN 

IDFEVAL  X/X 
END) 

1 1 .  Acknowledgement .  The  author  la  grataful  to  Prof.  Karl  Nickel  for  helpful 
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